Non-adiabatic dissociation of molecules and BEC loss due to shock-waves 
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Recent experiments have shown the likely appearance of coherent BEC atom-molecule oscillations 
in the vicinity of a Feshbach resonance. In addition, a new loss mechanism was observed, whereby 
the loss of atoms from the BEC is inversely dependent on the rate of change of the applied magnetic 
field. We present here a phenomenological model which gives a good description of the scaling 
properties of this new decay process, by attributing it to non-adiabatic dissociation of molecules by 
a propagating shock-wave. The model has only two free parameters, which specify the size of the 
" shocked-region" , and can be readily tested by future experiments. 
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In a recent series of papers the coherent creation of 
molecules from an atomic BEC close to the Feshbach 
resonance, was demonstrated jl], g|. Both sets of exper- 
iments have the following principle characteristics. In 
these experiments a stable 85 Rb BEC (at initial field B. L ) 
was subjected to a well controlled linear magnetic field 
pulse which brought it over a certain rise time t r to a fi- 
nal field Bf, close to the critical field Bq ~ 154. 9G of the 
Feshbach resonance. After a certain holding time th the 
field was ramped back to the initial value and the BEC 
was imaged. There was the distinct periodic appearance 
of a hot (^150nK) burst of atoms, and an unseen com- 
ponent, all of which remained coherent with the atomic 
BEC component. The unseen component is presumably 
the molecular component, since the Rabi frequency for 
the periodic change agrees with the energy difference 
between the free atoms and the bound molecular level 

Earn(Bf) @J. 

In addition to the periodic decrease and revival of the 
BEC, the results indicate an new mechanism of loss, 
which is seen to be much faster than the atomic two- 
and three-body inelastic loss rates All attempts to 
explain this loss in terms of adiabatic dynamics, treating 
the BEC as isotropic and uniform, have not been suc- 
cessful f|, ||- It is therefore useful to try a different ap- 
proach. The decrease in the atomic BEC could be due to 
molecules being formed or dissociated non-adiabatically, 
so that they do not take part anymore in the coher- 
ent oscillations jl). One possible non-adiabatic mecha- 
nism is the creation and propagation of shock waves in 
the atomic-molecular cloud. Shock waves in a classical 
gas are supersonic discontinuities in the temperature, en- 
tropy, pressure and density q], with all these quantities 
being increased by the passing shock-wave. Similar dy- 
namic discontinuities in the trapped cloud can result in 
non-adiabatic dissociation of molecules and therefore loss 
of coherence of the resulting atoms with respect to the 
atomic BEC. It is this process which we shall describe us- 
ing a phenomenological model here. We will not attempt 
a full simulation of the dynamics of the atomic and molec- 
ular clouds, but will describe the scaling of BEC loss due 
such a process. 



We begin by giving a qualitative description of mecha- 
nisms which can cause the appearance of shock-waves in 
the trapped cloud, due to the rapid change in the applied 
magnetic field. When the magnetic field is changed from 
its initial value Bi ~ 166G to its final value Bf close 
to the Feshbach resonance, the elastic atomic scattering 
length a changes by a factor of up to (2) ~ 4000, given 

by 

a = abg _ ~B~B~^) (1) 

with aft g = — 450ao, ao = 0.053nm and A = 11G. 

The equilibrium value of the density, using the usual 
Thomas-Fermi limit @, is: n eq oc a~ 3/5 . The BEC there- 
fore finds itself with the initial density n, which is much 
larger than the new equilibrium density n eq that corre- 
sponds to the new scattering length a(Bf). This means 
that the BEC cloud wants to greatly expand in response 
to the increase in a. The new equilibrium chemical po- 
tential /i eq = ATTh 2 n eq a/m oc a 2 / 5 is also very different 
from the chemical potential at which the cloud finds itself 
prior to any expansion, i.e. [i = 4irh 2 na/m. 

The expansion of the cloud to the new equilibrium 
size happens over a relatively long time scale, given by 
the largest trap frequency, which in the experiments is 
Ttrap — 30ms. But even before the cloud has time to 
readjust, the density begins to change through the propa- 
gation of density (sounds) waves, similar to the case of an 
expanding classical gas. Rarefaction waves at the BEC 
surface converge on the axis of symmetry and reflect as 
a shock- wave (Fig.l). Referring to a detailed simulation 
of a similar dramatic change in the scattering length || , 
we note that the imploding dynamics of the BEC create 
localized spikes of density (estimated to be up to lOn), 
of lateral size || W ~ 0.1/im. The width of the "shocked 
region" W is an essential parameter in our model (Fig.l). 

Another source of perturbation which can cause the 
shock-waves to appear, is the coherent conversion of 
atoms into molecules. This process will not be uniform 
in amplitude, being denser in the middle of the cloud, 
or in its frequency, due to the different local magnetic 
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field across the cloud. The BEC will therefore have a 
nonuniform density of molecules, oscillating (on average) 
at the Rabi frequency 0, |) 10 — 40kHz, depending on 
the holding field Bf. This time scale of 100 — 20/j,sec is 
of order of the time-scale of the observed loss rate. 

Our premise is therefore the following: when the mag- 
netic field is suddenly changed, acoustic waves propagate 
in the cloud in an attempt to restore the uniformity of 
the chemical potential, which is nonuniform due to the 
density profile of the cloud. Additionally, the BEC will 
be perturbed by an oscillating and nonuniform density of 
molecules. The consequent fluctuations (waves) in den- 
sity will result in reverberating shock-waves. Note that 
by the term shock-waves, we mean in general a propa- 
gating, non-adiabatic excitation. 

We estimate the amplitude of the shock-waves to be 
proportional to the gradients of chemical potential which 
form across the BEC due to the finite time it takes an 
acoustic wave to traverse the relevant region of the cloud. 
The discontinuity in the chemical potential at the shock 
front A/z, is therefore estimated to be the accumulated 
change in (i over the time r = W/v s it takes a sound 
wave to cross the " shocked-region" of size W, with the 
velocity of sound given by the standard expression [J7J: 
v s = fi/m (where we use the average density (n) and 
the velocity of sound at the holding field Bf). Note 
that weak shock-waves propagate at the speed of sound 
||. The chemical potential change, accumulated over the 
time t is 

A/Lt(r) = h(t) — /io — 47r/i 2 na(r)/m (2) 

(where fi = 4irh 2 nao/m, ag = (1 — A/(£>; — -Bo)) — 
— 2A), where 

a{T) = °>» (' - m^w) 
mr)-{l>; (3, 

describing the linear change in the magnetic field with 
time j^]. The dependence on the rise time t r is there- 
fore clear: the faster the rise, the larger are the chemical 
potential gradients created in the "acoustic" time r, and 
the resulting amplitude of the shock-waves. 

Note that our treatment is inherently dependent on 
the inhomogeneity of the trapped cloud (Fig.l); in an in- 
finitely uniform cloud there will be no shock waves. The 
process of shock-wave creation, even though we do not 
describe it here explicitly, depends on inhomogencitics 
which propagate in the cloud. We further immediately 
conclude that the decay time-scale r is largely indepen- 
dent on the BEC density n, as was found @. This behav- 
ior arises if the width of the " shocked-region" W and its 
length, are proportional to the equilibrium radius of the 
initial BEC, which is tx n 1 / 2 , while the sound velocity is 
also v s oc n 1 / 2 j7j. 

The sound velocity is given by the magnetic field dur- 
ing the hold time, i.e. closest to the Feshbach resonance. 



For Bf = 156. 7G we have v s — 8mm/sec, for initial BEC 
density (n) = 1.9x 10 13 cm~ 3 pL To compare our calcu- 
lation with the measured data [||, we used the measured 
decay time [||: r = W/v s — 13.2^sec, which corresponds 
to W ~ 0.1/xm. 

We now propose that these shock waves excite the 
molecules non-adiabatically and create pairs of atoms, 
which are therefore lost from the BEC (a similar mech- 
anism of molecular excitation/decay was mentioned in 
||). For the dissociation probability per molecule, due 
to the passing shock-wave, we assume a Boltzman-likc 
factor: p(r) — exp(— E am / A[i(t)). In a classical gas a 
passing shock front heats the medium, and the energy 
jump at the front, Afj,, plays the role of effective temper- 
ature |^| . Assuming that the shock- waves propagate over 
the width W at constant speed v Sl and cover a constant 
fraction / of the cloud (i.e. the "shocked-region"), the 
rate of loss of BEC is given by 

dNf = Nfpjr) 
dt t 

N(t r )/N = / e -W(r)/r + (!_/) ( 4 ) 

where we took the total time for the shock-wave deple- 
tion process to be empirically t to t = th + t r /4 d]> Nf is 
the number of atoms in the " shocked region" , which ini- 
tially is equal to fN . The parameter / is introduced in 
order to take into account the observation that the loss 
mechanism we describe here affects only a finite fraction 
/ ~ 90-60% of the cloud (Fig.2 of §). A possible reason 
for this behavior could be that the elongated shape of the 
cloud makes the conditions for effective shock-wave cre- 
ation and subsequent molecular dissociation appear only 
along its central part (Fig.l). This effect could explain 
why the proportion of the cloud which is affected by this 
loss mechanism is smaller in the case of the cloud with the 
smaller density [pi. This BEC will be shorter and there- 
fore have a proportionately shorter central " shocked" re- 
gion. The two free parameters of our model are therefore 
the width W (which determines r) and the relative num- 
ber of atoms / in the " shocked-region" , out of the total 
number of atoms in the cloud. 

To take into account the background processes of two- 
and three-body inelastic loss rates we next multiply 
N(t r )/N (§) by: exp [-t h /T 2 - 3 ], with a typical decay 
time of T2_3 ~ 0.12 x 10~ 2 scc (in the vicinity of the 
Feshbach resonance Q). 

In Fig.2 we compare our calculated fraction of remain- 
ing BEC N(t r )/N (Eq.(|)) as a function of t r , for dif- 
ferent hold times th, with the experimental data In 
the limit of vanishing rise time t r , the fraction of the re- 
maining BEC shows an exponential decay with the hold 
time th (Fig.2 of Q). For increasing rise time we first see 
the increasing loss, as long as t r < t, simply due to the 
increase in the overall time spent close to the Feshbach 
resonance. As soon as t r > t, the loss is decreased, due 
to weakening shock- wave amplitude (^|^). The overall 
agreement is very good. The discrepancy at the longest 
holding time (th = lOO/isec) arises from the fact that 
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the reverberating shock- waves weaken with time, making 
longer holding times relatively less effective. This effect 
is not taken into account in Eq.(^), where we assumed 
that the dissociation probability p(r) is constant with 
time. We are therefore overestimating the loss of BEC 
for the longer holding times. A much better agreement 
is achieved for the case of th = 100/isec, if an effectively 
reduced value of th ~ 50/isec is used in the calculation 
, describing the decay of the shock-waves (dashed-line 
in Fig. 2). The background decay due to two- and three- 
body processes |^] is shown by the dotted line. 

In Fig. 3 we compare our calculation (Q), for a constant 
holding time th = l^sec (Fig. 4 of ||), and varying final 
field Bf. The velocity of sound, and therefore the decay 
time t, both depend on the magnetic field Bf, with r 
decreasing as we approach the resonance field Bq . Addi- 
tionally, the dependence of the dissociation energy E am 
on the magnetic field was previously measured [Q . Again 
we find a good agreement, except for the longer rise times 
and closest approach to the resonance. At this field, our 
empirical assumption of t r /A being added to the over- 
all decay period is questionable, with a larger proportion 
probably active. On the other hand, this means that the 
shock-waves have decayed for a longer time too, resulting 



in our overestimation of the BEC loss, as discussed above 
in relation to Fig. 2. 

To conclude, we have presented a phenomcnological 
model which attributes the loss of atomic BEC to non- 
adiabatic dissociation of molecules by imploding shock- 
waves, created as a result of the rapid change in the mag- 
netic field. This simple model appears to capture the 
main physical mechanism at work, as it describes the cor- 
rect dependence of the observed BEC loss on the various 
physical parameters. Rigorous numerical simulations are 
needed in order to substantiate this proposal. It could be 
interesting to test this model further by repeating the ex- 
periments with BEC clouds of different geometries (e.g. 
prolate vs. oblate), and by producing shock- waves in con- 
trolled regions of the cloud by applying time-dependent 
and spatially non-uniform electro-magnetic fields. 
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FIG. 1: A schematic illustration of our proposed scenario, 
whereby (a) imploding acoustic (rarefaction) waves are (b) 
converging on the axis of symmetry to produce a cen- 
tral " shocked-region" (dark region), where the non-adiabatic 
molecular dissociation takes place. 
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FIG. 2: The calculated (solid lines, Eq.(|)) and measured 
[Q] (symbols) remaining fraction of atomic BEC as a func- 
tion of the rise time t r , for different holding times th = 
1, 5, 15, 35, 100/^sec (from top to bottom). The dotted line 
gives the background decay due to two- and three-body pro- 
cesses The dashed line describes the calculation for a 
holding time of th = 50/isec. 
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FIG. 3: The calculated (solid lines, Eq.(|)) and measured 
M (symbols) remaining fraction of atomic BEC as a func- 
tion of the rise time t r , for different holding magnetic field 
Bf = 158, 157.2, 156.7, 156G (from top to bottom). The dot- 
ted line gives the background decay due to two- and three- 
body processes M. 



